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We study the entanglement properties of deconfined quantum critical points. We show not only that these 
critical points may be distinguished by their entanglement structure but also that they are in general more highly 
entangled that conventional critical points. We primarily focus on computations of the entanglement entropy of 
deconfined critical points in 2 + 1 dimensions, drawing connections to topological entanglement entropy and a 
recent conjecture on the monotonicity under RG flow of universal terms in the entanglement entropy. We also 
consider in some detail a variety of issues surrounding the extraction of universal terms in the entanglement 
entropy. Finally, we compare some of our results to recent numerical simulations. 



I. INTRODUCTION 

Conventional ground states of bosonic quantum matter are 
distinguished by the presence or lack thereof of broken sym- 
metry and associated long range order (LRO). These ground 
states have prototypical wavefunctions that are direct prod- 
ucts of local degrees of freedom, and hence contain only short 
range quantum entanglement. In the last several years much 
attention has been lavished on ordered states which do not fit 
the broken symmetry/long range order paradigm. The best 
studied and simplest examples are gapped topological phases 
with a low energy effective description in terms of a topolog- 
ical field theory as introduced in Refs. 1,2. More complex 
examples also exist in the form of gapless liquid-like phases 
which also do not fit the broken symmetry paradigm. A com- 
mon feature of all these 'exotic' states is that their prototypi- 
cal ground state wavefunctions are not direct products of any 
local degrees of freedom. These states have long range en- 
tanglement (LRE) in their ground state wavefunction. Thus 
the study of the entanglement properties of the ground state of 
various quantum phases and their associated phase transitions 
has become a topic of fundamental importance. Entanglement 
based ideas have led to a classification of gapped quantum 
phases in one dimension described in Refs. 3,4, efficient ways 
to characterize critical theories in one dimension described in 
Refs. 5,6, a partial characterization of topological phases in 
two dimensions as shown in Refs. 7,8, and much more. 

Due to these and other considerations, there has been a 
recent profitable exchange of ideas between quantum infor- 
mation science and quantum many-body physics. Broadly 
speaking, the study of entanglement entropy has provided us 
with a rough guide to the landscape of many-body states. It 
guides us to the physical region of Hilbert space where many- 
body states live. Most importantly, entanglement entropy is at 
present one of the best available quantities to concretely probe 
the structure of long range entanglement in quantum many- 
body ground states. In this paper, as part of the broader effort 
to understand the entanglement properties of zero tempera- 
ture phases and phase transitions, we describe the structure of 
entanglement at certain deconfined quantum critical points in 
higher dimensions with a special focus on two dimensions. In 



particular, we will compute universal subleading corrections 
to the entanglement entropy for a variety of deconfined quan- 
tum critical points in two dimensions described in Refs. 9,10. 

Deconfined critical points lie outside the scope of the Lan- 
dau paradigm for phase transitions because they are most nat- 
urally described in terms of "deconfined" degrees of freedom 
that exist at the critical point (but perhaps not in either phase). 
The emergence of these deconfined degrees of freedom and 
associated gauge fields suggests that the ground state of these 
critical points may have "stronger" long range entanglement 
as compared to conventional quantum critical points (i.e those 
described within the Landau paradigm). It is therefore inter- 
esting to explore their entanglement properties in some detail. 
Deconfined critical points are useful paradigmatic examples 
of exotic gapless systems with strongly non-quasiparticle like 
excitation spectra. Studying their entanglement structure is a 
step in exploring the relationship between the breakdown of 
the concept of a quasiparticle to describe the excitations of a 
many-body system and possible long range entanglement in 
the corresponding ground state wavefunction. 

A simple example of a deconfined critical point is pro- 
vided by a quantum phase transtion between a superfiuid and 
a Zi topological Mott insulator of bosons at integer or half- 
integer filling on a lattice (see Refs. 1 1-16). The latter phase 
has gapped fractional charge- 1/2 bosons (dubbed chargons) 
and gapped Z2 vortices (dubbed visons) which are mutual 
semions. The transition to the superfiuid goes hand in hand 
with the destruction of topological order. Despite the pres- 
ence of a natural Landau order parameter for the superfiuid 
phase, the critical theory is not described by a conventional 
Landau-Ginzburg- Wilson action in terms of this order param- 
eter. Rather it is most simply described as a condensation tran- 
sition of the fractionally charged boson in the D = 2 + 1 XY 
universality class. The exponents v and z are the same as in 
the usual XY universality class, but the anomalous dimension 
of the superfiuid order parameter is large (effectively because 
the superfiuid order parameter is a composite field). The uni- 
versality class of this transition has been denoted XY* to dis- 
tinguish it from XY. Coming from the superfiuid side the de- 
confinement of the fractional boson occurs already at the crit- 
ical point and continues into the Mott phase. This deconfined 
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critical point provides an excellent setting for us to contrast 
the entanglement behavior with conventional quantum critical 
points by simply comparing the universal subleading behavior 
of the entanglement entropy of the XY* and XY transitions. 
Indeed, we will show that this universal part is different for the 
two transitions with the difference precisely due to the emer- 
gent topological structure of the proximate Mott insulator. 

Another well studied example of a deconfined critical point 
is that of a continuous Landau-forbidden quantum phase tran- 
sition between Neel and Valence Bond Solid (VBS) phases 
of spin- 1/2 quantum antiferromagnets on a square lattice (see 
Refs. 9,10). The critical field theory is naturally described in 
terms of emergent gapless bosons that carry fractional spin 
and are coupled to an emergent gapless U(l) gauge field. 
In this case we argue that the universal part should be cor- 
rectly given by calculations in the non-compact CP 1 confor- 
mal field theory that describes the low energy physics - mi- 
croscopic lattice models differ from this field theory only by 
irrelevant operators which do not affect the universal contri- 
bution to the entanglement entropy. 

Further interest in these particular models comes from pio- 
neering numerical work reported in Refs. 17,18. They studied 
the entanglement properties of a simple realization of an XY* 
transition. We will analyze their calculations in the context of 
our results later in the paper. Even more recently, Ref. 19 
studied various universal terms in the entanglement entropy 
of the square lattice Heisenberg model. That work raised cer- 
tain puzzles associated, for example, with corners. Our work 
here suggests one possible resolution of those puzzles. They 
compare their Monte Carlo and series results for corners with 
that of two free bosons (a non-compact spin wave descrip- 
tion). However, it may be important to remember the com- 
pactness of the actual order parameter. In the case of an XY 
ordered phase, the compactness of the order parameter may 
be approached by coupling a non-compact free boson to a Z 
gauge field. As we will argue repeatedly in this work, such 
couplings typically have the effect of shifting universal terms. 
Our work below leads us to believe that the gauge field cou- 
pling would make the corner contribution more negative in 
closer agreement with the numerical results. At the very least, 
there could be a crossover length scale that would be have to 
be exceeded to truly detect the non-compact boson corner con- 
tribution. Simply speaking, we are reducing the number of de- 
grees of freedom and hence should reduce the entanglement. 
We also note that the J — Q model introduced by Sandvik in 
Ref. 20 seems to realize the deconfined Neel to VBS transi- 
tion, and since this model is sign problem free, it is also a very 
attractive target for investigations of entanglement entropy at 
deconfined critical points. 

Now we outline the details of our calculations. The en- 
tanglement entropy of a region R is defined as S(R) = 
—ten(PR m Pr) where pn is the density matrix for R. When 
considering multiple regions, we simply concatenate the la- 
bels, hence S{R\R2) is the entropy of the combined region 
Ri Ui?2- The basic rule governing the entanglement entropy 
of a local many-body system is the boundary law. In d spatial 
dimensions the entanglement entropy of a region R of linear 
size L is roughly L d_1 i.e. it is proportional to the size of 



the boundary of the region as reviewed in Ref. 21. Excep- 
tions to this rule include gapless systems in one dimension 
as shown in Refs. 22,23 and systems with a Fermi surface 
in higher dimensions as argued in Refs. 24-27, but the rule 
is obeyed by a wide variety of gapped and gapless phases 
in higher dimensions. It is a super-universal feature of lo- 
cal many-body systems. That being said, the prefactor of the 
boundary law term in the entanglement entropy is not directly 
meaningful because it contains details of the microscopic reg- 
ulator. Different microscopic realizations of the same phase 
or phase transition can have different values for this prefactor. 
In other words, while the existence of the boundary law term 
is super-universal, its precise value is not universal at all. 

All this information can be conveniently summarized by an 
intuitive scaling argument described, for example, in Ref. 28. 
To obtain the entanglement entropy of a region R we take a 
renormalization group perspective and add up the contribu- 
tions to the entanglement entropy from each scale. The key 
assumptions are locality in space and approximate locality in 
scale (so that we may add contributions from different scales). 
Let r indicate the length scale of interest in our many-body 
system, r runs from the ultraviolet cutoff e (the lattice scale) 
to oo in the long wavelength limit. Spatial locality suggests 
that the contribution to the entanglement entropy at scale r is 
proportional to the size of the boundary measured in units of 
r: dS(r) ~ (L/r) l . We integrate this contribution from 
the cutoff e to the minimum of the the correlation length £ 
(gapped) and the region size L (gapless). The appropriate 
renormalization group measure is d\nr = dr/r. The result 
of this simple calculation is a non-universal boundary law for 
gapped phases in d = 1, a universal logarithmic dependence 
for gapless phases in d = 1, and a non-universal boundary 
law for gapped and gapless phases in d > 1. Thus as we saw 
above, to extract universal features in higher dimensions we 
must turn to sub-leading terms in the entanglement entropy. 

The general structure of sub-leading terms can be un- 
derstood via a generalization of the scaling argument given 
above. In Ref. 28 a "twist field" formalism was used to inves- 
tigate the short distance structure of entanglement. The result 
of this expansion is that only even or odd terms in an L/r 
expansion appear as appropriate; furthermore, the pattern ob- 
tained in Ref. 28 agrees with known results (see Ref. 29,30) 
in 1 > d < 3 spatial dimensions. In two dimensions, the en- 
tanglement entropy of a smooth region of linear size L in a 
gapped phase has the form S(L) = aL — /3 + .... The co- 
efficient a is non-universal, however, the coefficient (3 is uni- 
versal and partially characterizes the topological order of the 
phase as shown in Refs. 7,8. [3 = for any phase adiabatically 
connected to a trivial product state e.g. a bosonic Mott insula- 
tor at integer filling. For topological phases characterized by 
a modular S-matrix Sg, we have j3 — \iiSq = — ln£> with T> 
the total quantum dimension, a result we will re-derive later. 
In a certain precise sense the coefficient f3 increases with the 
quantum complexity of the topological phase. Gapless phases 
in two dimensions have a similar form for their entanglement 
entropy, but other terms are also possible, for example, cor- 
ners will give logarithmic terms in the entanglement entropy 
(see Ref. 31 for a sample calculation and Ref. 28 for a scaling 
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argument). Note that there is some subtlety in giving meaning 
to f3 in the presence of such logarithmic corner terms. Fur- 
thermore, in a gapless phase (3 will depend on region shape, a 
fact we will sometimes emphasize by writing j3 = (3(R). For 
example, a long thin strip will have (3(L,w) ~ L/w with L 
the length and w the width of the strip. It is thus important to 
give a precise protocol for defining and extracting (3. We will 
return to these questions below. 

Before moving on, we wish to comment on a slightly 
counter-intuitive feature of the above discussion. We claimed 
that phases with larger (3 are more entangled, but looking at 
the formulae above, it seems that a larger f3 corresponds to 
a smaller entanglement entropy. We may understand this as 
follows. For any large smooth region in a gapped phase, the 
entropy has the form S = aL — j3. We may change the value 
of a by adiabatically deforming the Hamiltonian within the 
phase, but the value of f3 is fixed and universal. Thus because 
/3 contributes negatively to a positive quantity, the entangle- 
ment entropy, we know that in any phase with non-zero /?, 
we cannot reduce a to zero by deforming the Hamiltonian. 
In other words, a phase with non-zero j3 has no microscopic 
Hamiltonian in which the ground state is a product state. As (3 
increases, the "obstruction" to a product state also increases, 
thus establishing that larger f3 corresponds to a "more entan- 
gled" phase of matter. 

Our basic result for the XY* transition is that for simply 
connected regions (3xy* = Pxy + Pz 2 where /3z 2 = In 2 
is the topological entanglement entropy of a Z 2 gauge theory. 
As described above this simple result explicitly shows how the 
universal entanglement entropy acquires an additional topo- 
logical contribution not present at the XY critical point. This 
gives precise meaning to the notion that this deconfined crit- 
ical point has stronger entanglement than the corresponding 
conventional critical point. In addition to the simple example 
of the XY* transition, we obtain similar results for an infi- 
nite class of similar deconfined critical points. These decon- 
fined critical points consist of bosons undergoing a superfluid- 
insulator transition where the insulator is an exotic fractional- 
ized phase. The description is in terms of emergent bosonic 
degrees charged under deconfined discrete gauge fields. We 
also argue that the entanglement entropy of lattice models re- 
alizing the deconfined critical point separating Neel and VBS 
phases has the same universal properties as the critical CP 1 
model. We also explain why the conjecture of Ref. 32 con- 
cerning the existence of an entanglement c-theorem shows 
that the CP 1 model is more highly entangled than the usual 
0(3) fixed point. 

This paper is organized as follows. First, we describe the 
XY* transition and compute the entanglement entropy di- 
rectly from a model wavefunction. Second, we generalize this 
argument to a large class of models using a partition function 
argument recently described in Ref. 33. Third, we discuss 
the properties of the CP 1 model and its field theory. Fourth, 
we discuss a variety of procedures for defining and extract- 
ing f3. Finally, we conclude with a discussion of the broader 
relevance of our results. 



H. XY* TRANSITION 
A. Description of transition 

Consider a system of bosons on a lattice at a commensurate 
filling of, say 1/2, per site on average in two space dimensions 
described by a Hamiltonian with the general structure 

" XXVi^ ■ !>■<■) ■ - (i) 

with ip\ being the creation operator for a boson at site i. 
For concreteness we will specialize to a square lattice though 
our results apply equally to other lattices. If the boson hop- 
ping dominates over all interaction terms a uniform superfluid 
phase will result. With increasing strength of short ranged 
repulsive interactions a variety of Mott insulating phases re- 
sult. Some simple possibilities are Mott insulators that break 
lattice symmetry such as checkerboard density wave or bond- 
centered stripe insulators. Our focus will however be on Mott 
insulators that preserve all symmetries of the Hamiltonian. 
For a square lattice and at half-filling such Mott insulators 
are necessarily exotic (see Refs. 34,35): they either have 
topological order and a gap to all excitations in the bulk or 
some even more exotic order and gapless excitations. We 
will mainly restrict attention to the simplest such possibility, 
namely a gapped Mott insulator with topological order that 
is described by the deconfined phase of a Z 2 gauge theory. 
Such a phase has gapped fractionalized excitations (dubbed 
chargons) that carry 1/2 the boson charge. In addition here 
are gapped charge-neutral vison excitations that are mutual 
semions with the chargons. For various examples of micro- 
scopic boson models that realize such a phase and the associ- 
ated phase transitions, see Refs. 12,13,15,16,18,36. 

Being topologically ordered the fractionalized Mott insu- 
lator has Long Range Entanglement in its ground state that is 
partially captured by the concept of Topological Entanglement 
Entropy. Recent numerical work has successfully demon- 
strated this property for a specific model. Here we are inter- 
ested in the quantum phase transition between the superfluid 
state and the fractionalized Mott insulator as the interaction 
strength is varied at fixed filling. This transition is conve- 
niently described as a condensation of the chargons while the 
vison gap stays non-zero. Specifically if we describe the char- 
gons by a field b, the transition is described by the action of 
the space-time D = 3 dimensional XY fixed point Sxy[b]- 
Despite this however the transition is not quite in the 3D XY 
universality class. This is because the field b itself carries Z 2 
gauge charge and hence is not directly observable. The phys- 
ical boson operator tp ~ b 2 and hence has a large anomalous 
dimension r} « 1.4. Indeed the only allowed physical oper- 
ators at the fixed point are the subset of the operators at the 
3D XY fixed point which are invariant under the local Z2 
gauge transformation associated with changing the sign of bi 
at any given site. The fixed point describing the superfluid- 
fractionalized Mott insulator transition is thus distinct from 
the 3D XY fixed point and is denoted as an XY * fixed point. 
This 1 1-16 and other similar transitions 37-40 have been de- 
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scribed in a number of contexts in previous papers. Most re- 
cently this transition has been explored numerically in Ref. 
16 and the predicted large rj and other properties have been 
demonstrated. 

As the simplest example of a non-Landau transition out of 
an ordered phase the XY* transition provides a good opportu- 
nity to examine the entanglement structure of an exotic phase 
transition. Our goal will be to compare the universal structure 
of the long range entanglement at the XY* fixed point with 
that of the XY fixed point. This task is simplified by noting 
that as far as the universal properties of the transition are con- 
cerned we might as well take the vison gap to infinity. As the 
visons correspond to the magnetic flux of the gauge field this 
means that the 'electric' field lines of the gauge field can be 
taken to be tensionless. 



B. Exact wavefunction computation 

Consider an L x L square lattice with bosons bi on the sites 
and Z 2 gauge fields <r|j- on the links. The bosons are charged 
under the gauge field and we take the gauge dynamics at the 
lattice scale to be tensionless. The Hamiltonian is 



h.c. — /i > n,: 



-U 



<i]> 

ni(ni - 1) 



K 



E 

ijkl£d 



(2) 



with i%i = b^bi. We also impose the constraints Ci 



ii J rx G ii—x u ii+y L> ii — y 



1. To study the XY* transition 
we tune the boson filling to be one per site and study the super- 
fluid to insulator transition in the presence of the gauge field. 
Because we are in the tensionless limit, the gauge fields cr~ 
are constants of the motion. Choose a gauge where crL = 1 
for all ij and call the ground state of the boson model \B). Im- 
posing the constraints attaches the ends of electric field lines 
to the bosons. 

Write the boson state as \B) = ~^2 p c p \B p ) where p is an 
index that records the even/odd parity of n, at each lattice 
site. Thus p runs over 2 2i " values. We take (B v \B q ) = 
Spg, The index p records e l7!,li for each site and hence 

Pu+xOu-xOfi+yOfi-y for each site. There are 2 2L " states 



in the unconstrained gauge space. For each configuration of 

(^ii+x^ii-x^ii+y^ii-y there are 2 L * states. These states can 
be described as the set of all string states with strings ending 
on the sites specified by p. They can be constructed by starting 
with a particular state with the right string endings and acting 
with plaquette operators to generate all other states. Each state 
in the sum has an equal weight because of the tensionless na- 
ture of the strings. 

Imposing the constraint is simple in the p basis. We simply 
replace \B)\G) with 



Y,c P \B p )\G p ) 



where \G P ) are the states described above. 

Now we must ask about the spatially partitioned entangle- 
ment entropy. In fact, we will focus on the second Renyi en- 
tropy. Partition the system into two regions i?i and R2 in the 
L — >• 00 limit. To be concrete we can take region R\ to be a 
disk and region R2 is the rest of the universe. We will compute 
S 2 (Ri). 



1. Bosons only 
For each \B p ) we make a decomposition 

Xp 

\B p ) = £ ^a p \B pap )l\B pap ) 2 
a p = l 

Just to be clear, we have a different set of A's for every p, 
but in the interest of avoiding a notation explosion, we try to 
suppress this fact. The subscript on a p indicates the value of 
p. We may take the A's to be real. The full density matrix is 
pure p — \B)(B\. The reduced density matrix for region 1 is 



Pi 



-^c p c* £ ^a p ^a q \B pap )(B qaq \ 1 (B pap \B qaq ) 2 



I'll 



In general, the only thing we can say about (B poip \B qoiq ) 2 is 
that it contains a factor of S P2q2 where p 2 and q 2 indicate the 
parts of p and q specifying information in R 2 . 



We then need \ii(p\) to compute S 2 {R\). We find 



tn(p?) = £ 



pqrs 



^ol p Aa q ^a r Aa s \^pa p \^qa q /2 \^ra. p \ ^sa q ) 2 {^qa q \ Brat r )l (^sa s \ ^pa p ) 1 



(3) 



a p a a a r a s 



Again, all we can say about this expression is that (B pctp \B qa ) 2 contains a factor of 6 P2q2 and (B pap \B qa )i contains a factor of 
8 Piqi where p and q could be any of p, q, r, s. Note the usual replica branched surface kind of structure present in this otherwise 
messy looking formula: p goes with q and r goes with s in R 2 while q goes with r and s goes with p in R\. 



2. Bosons + gauge fields 

We now need a similar decomposition for the gauge field 
states. Write 



2l°«i 



\G V 



E 



i 



\G p i)i\G p i) 2 



There are important simplifications in this case. First, the 
Schmidt rank is the same for all p. It is always 2^ dRl ^ 1 be- 
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FIG. 1: Schematic of the computation of tr(p?) illustrating the dif- 
ferent connectivities in Ri and i?2- Vertical bonds represent contrac- 
tions of physical degrees of freedom. 



cause we may specify the presence or absence of a string on 



each link piercing the boundary except for an overall global 
constraint requiring the number of strings entering to be even 
or odd depending on the number of charges inside R\ (as de- 
termined by p). Similarly, we may imagine for every p that 
A just labels precisely these boundary string configurations. 
Thus given two choices p and q that agree on R2 (but may 
disagree on i?i), the states from R2 appearing in the decom- 
position above are actually the same. In other words, we have 
(G p i\G q j)2 = $ijS P2 q 2 . This is a special property of the 
gauge fields in this theory that permits progress to be made. 

The full boson/gauge reduced density matrix is thus 

Pi=5I c p c 9 ^a P ^a q \B pap ){B qaq \i{B pa \B qaq ) 2 {4-) 

Pq OlpOtqU 

x '^jaRTpT I G p i)(G qJ 1 1 (G pI I G qJ ) 2 



Now again we compute tri (p\) 

= C P C q C r C *s K p K q K r K s (B pap \B )2(B roiv \ B solq )2(B qaq \B rar )i(B SOCs \ B pap )i 



(5) 



pqrs 



a. v a. Q ct r a. s I J K L 



1 



(G p i\G q j)2(G r K\G sL } 2 (G q j\G r K)i(G S L\G p i) 



2l9«i|-i ; 

Now we use the fact that (G p i\G q j)2 = &i,j5 P2q2 and similarly for inner products in R\. 

ttlCPl) = C P C q C r c * ^a p ^u q ^a r ^aAB p a p \B qaq )2(B ra \B saq )2{B qaq \B rar )i{B sa \B pap )i 

2 



(6) 



pqrs 



a v a Q a. r a s I J K L 



1 



SljS p2 g 2 5KL$r 2 S2 SjK&qin 8li5s iPi 



The sums over IJKL now reduce to a sum over I = J = K = L giving a factor of 2^ 9Rl ' x . Moreover, the explicit delta 
functions may be absorbed into the corresponding boson overlaps e.g. {B pap \B qaq )28 P2q2 = (B pap \B qaq )2 since (5 P2q2 ) 2 = 
5 P2 q 2 and (B pap \B qa )2 already contains such a factor. Thus we find 



tr l(Pl) — 5^ C p C*C r C* ^a p ^a q ^a 7 .^a s (B pap \B qaq } 2 (B rap \B solq } 2 {B qaq \B rar )i(B sas \B pap )i 



a p a a a r a s 



1 



2l3fli|-i 



(7) 



osons only 



This final line implies that 5*2 ) = ^(^Obosons only + m S assoc i ate d with a critical point), hence there is a precise 
(|flRi I - 1) In 2. In fact, this relation holds for all Renyi distinction between these two contributions to p XY . . 
entropies and hence holds for the entanglement entropy as 
well. Note that in this case the gauge field contribution is 
in fact independent of n. Thus if we parameterize the en- 
tanglement entropy at the critical point as S — aL — f3 
then we have Sxy = a L — (3xy, S% 2 = a> ' L — /?z 2 , aR d 



C. More complicated regions 



Sxy* = ot"L - I3 X y* with j3 X v = Pxy + , 
that (3z 2 does not depend on region shape while @xy will (be- 



We have just seen via an explicit wavefunction computa- 
We note tion that at the critical point the entanglement entropy has the 



form aL — j3 with f3 = fixY + PZ2. However, we only es- 
tablished this result for a simply connected region. The alter- 



FIG. 2: Four regions defining the Levin-Wen procedure for extract- 
ing the topological entanglement entropy. The regions are chosen so 
that all boundary and corner terms will cancel. 



native derivation to be presented below based on a partition 
function argument also only tells us the answer for a simple 
disk geometry. The question thus remains, how does (3 depend 
on region shape and topology? As we will argue in a later sec- 
tion, the gapless contribution to j3 has both non-trivial shape 
dependence as well as non-trivial topological dependence. On 
the other hand, the topological contribution, being associated 
with local constraints, depends only on the topology of the 
region. 

In a gapped phase of matter in two dimensions, the entan- 
glement entropy of a region R of linear size L takes the form 
S = aL — (3(R) = aL — f3(disk)Ng where Ng is the num- 
ber of disconnected boundaries of R. This formula is valid up 
to exponentially small corrections in L/£. Corners will not 
contribute logarithmically growing contributions for L ^> £, 
however they can contribute terms of the In £/e which must be 
separated out from the smooth part /3. If we can arrange, in nu- 
merical computations, say, to study regions without corners, 
then j3 can be extracted from finite size scaling of S(L)/L. 
If corners cannot be avoided, as is often the case on a lattice, 
then extra care must be taken to remove the corner contribu- 
tions. 

Protocols to achieve this separation for gapped phases were 
presented in Refs. 7,8. In the case of Ref. 8, the entangle- 
ment entropy of four regions are subtracted in a certain way 
to remove all boundary law and corner terms. The four re- 
gions, shown in Fig. 2, call them Ri (i = 1, 4), look like 
an annulus (i?i), two half moons (R2 and R3), and two dis- 
connected strips R4. All regions are topologically collections 
of disks except for the annulus R±. The prescription is then 
(S(R 1 )-S(R 2 ))-(S(R 3 )-S(R 4 )) = (-2/3 + /3)-(-/3 + 
2/3) = -2(3. 

We may also ask, in the context of Z 2 gauge theory, say, 
what happens when either the electric or magnetic charges 



FIG. 3: A smoothed version of the annulus entering the Levin- Wen 
protocol. We need the structure of the reduced density matrix of this 
geometry as a function of the radii a and b. 



are allowed to fluctuate. This question may be given a sim- 
ple answer in the context of the toric code by putting one of 
the constraint terms to zero. There is now a macroscopically 
degenerate space of ground states, but we may still extract a 
"topological entropy" via the same procedure (volume, area, 
and corner terms are all subtracted out). The important point 
here is that the constant term no longer depends on the num- 
ber of disconnected components in a simple way. For exam- 
ple, two disconnected strips still give —2/3 = —2 In 2, but the 
annulus now only gives —j3 = — In 2. The prescription now 
gives — /3, half the original result. The interpretation of this 
result is that half the constraint structure has been destroyed 
with the other half remaining intact. 

Now we ask what happens in a gapless state. In the case 
of interest, electric charges are allowed to fluctuate, but their 
motion is not completely washed out (no large number of de- 
generate states). First, (3 will now depend on region shape in a 
non-trivial way. For example, a long thin strip (with smoothed 
corners) of length L and width w would have (3 ~ L/w. Also, 
corners will now give a logarithmic contribution of the form 
f(6) In (L/e) with 8 measuring the deviation from a straight 
line and / an even function. We also note that these consid- 
erations apply for regions embedded into an otherwise infi- 
nite system. Finite size corrections will modify all these state- 
ments, and some considerations along these lines will be rele- 
vant when comparing to numerical data. Already the prescrip- 
tions in Refs. 7,8 are less appealing because, although they 
subtract the corner dependence in a sensible way, the shape 
dependent parts of (3 are no longer guaranteed to nicely can- 
cel. We return to this point later. 

For now, we want to argue that the topological part of (3 for 
the annulus is 2 In 2 even at the critical point. Such a result 
is quite different from the answer obtained by putting one of 
the constraint terms in the toric code to zero. Since we are 



7 



no longer dealing with a gapped phase, we must specify the 
geometry of the annulus more carefully. Let a be the inner ra- 
dius and b the outer radius. We consider a limit of a, b — >• oo 
with a/b — > e.g. the outer radius is going to infinity faster 
than the inner radius. Call A the inner disk, B the annulus, 
and C the rest of the system as shown in Fig. 3. Because 
the state of the whole system ABC is pure, the entanglement 
entropy S(B) of the annulus is equal to S(AC). The den- 
sity matrix for this composite system is expected to factorize 
in the limit we consider giving S(AC) S(A) + S(C). To 
quantify this expectation, we turn to the mutual information 
I(A,C) = S(A) + S(C) - S(AC). The mutual informa- 
tion bounds connected correlations according to (OaOc)c < 
\\Oa\\ 2 \\O c \\ 2 I(A,C) as shown in Ref. 41 , and although we 
do not prove it, we expect this bound to be tight at the critical 
point. 1 Thus the mutual information decays as fast as square 
of the slowest decaying connected correlation function of nor- 
malized gauge-invariant operators in A and C. This means the 
density matrix pac factorizes as a, b — > oo as a power law in 
the separation b — a. Using this information, we can conclude 
that the shape dependent part of /? is simply that of two disks, 
up to power law corrections. This also suggests that the topo- 
logical parts is also equivalent to that of two disks, but given 
the importance of this claim, we further justify it. 

Now we turn to the topological part. Let us imagine fix- 
ing, as in the wavefunction calculation above, the distribu- 
tion of Z2 charges inside the annulus. Suppose there are an 
even number of charges. Without gapless charges, we would 
be forced to say that an even number of electric strings en- 
ter from both A and C, a situation we call (0, 0). However, 
with fluctuating charges we could also have an odd number 
of charges in A and C. Now there will be an odd number of 
strings entering from both A and C, but still an even number 
of strings entering the annulus through all its boundaries taken 
together. This situation we call (1, 1). In the toric code with 
constraint removed, the existence of these two sectors reduces 
f3 from 2 In 2 to In 2. The important question is, does the state 
in region B care if we have an even or odd number of parti- 
cles in A and C. Is the boson wavefunction (before the strings 
are attached) more like (i) (|(0, 0))ac + |(L l))Ac)l even )B 
orf//j(|(0,0)) A c|even,(0,0)) s + |(l,l))Ac|even,(l,l)) B ). 
Case ( i) is typical of the superfluid phase where bosons may 
fluctuate freely. In such a phase the the correlator (b^bj), 
containing an operator that moves us from (0, 0) to (1, 1), is 
finite because the state is long range ordered. This correla- 
tor is also finite in the cartoon state (i). Case (ii) is typical 
of the Mott insulating phase where the bosons are locked in 
place. In such a phase the correlator (6,4 is zero, as it is in 
the cartoon state ( ii) (assuming the states in B are orthogonal). 
Because at the critical point the correlator ipA be) decays with 
a universal power law, we argue that the critical point is more 
like the Mott insulator than the superfluid phase. Hence we 



One situation where the mutual information does not decay appears when 
we consider "cat states" in symmetry breaking phases. Such states can be 
relevant when considering finite size numerical simulations. 



should find that the topological part of j3 is 2 In 2. We empha- 
size again that this is all true up to power law corrections at 
the critical point. 

We briefly comment on the situation at finite temperature to 
illustrate the scales involved. We imagine starting at T = in 
the Mott insulating phase, then moving up to encounter first 
the quantum critical region, then a higher temperature region 
where the boson charge dynamics are disordered, and finally 
to an even higher temperature region beyond vortex gap. The 
topological entropy, defined by the usual protocol, should be 
—2 In 2 in the insulating phase and in the quantum critical re- 
gion as we just argued. Above the dynamical scale set by 
the boson motion, we should find — In 2 as is appropriate for 
the electrically disordered toric code. Finally, the second con- 
straint is lifted above the vison gap, and we find no topological 
entropy. 

III. PARTITION FUNCTION ARGUMENT 
A. Conformal transformation 

Now we give a more general argument for the structure of 
the entanglement entropy by mapping the entanglement en- 
tropy of a disk onto a problem of computing the partition 
function of a sphere, see Ref. 42. This transformation works 
whenever the theory in question is a conformal field theory. 
This includes topological field theories (a highly degenerate 
case of conformal field theory) and most z = 1 scale invari- 
ant theories. However, we should be careful applying these 
results to z 7^ 1 scale invariant theories. 

We want to compute the entanglement entropy of a disk 
B 2 in d = 2 dimensions, or more generally, a ball B d of 
radius R in d spatial dimensions. We begin by writing d + 1 
dimensional Minkowski space metric in spherical coordinates 
as 

ds 2 = -di 2 + dr 2 + r 2 d^_i, (8) 

where dVt 2 d l is the unit round metric on the d — 1 sphere. In 
the case of d = 2, we have d£l 2 ,_ 1 = d<f> 2 with cf> ^ <f> + 2ir. 
Now consider the coordinate transformation given by 

cos 9 sinh t/R 
1 + cos 9 cosh t/R 



1 + cos 9 cosh t/R 
The metric in these new coordinate now reads 

ds 2 = F 2 (— cos 2 9dT 2 + R 2 (d9 2 + sin 2 9dn 2 d _ 1 )) 

F = (1 + cos 9 cosh t/R)' 1 . (10) 

The conformal factor F 2 may now be removed by a conformal 
transformation. 

In these new coordinates the "entanglement Hamiltonian" 
is nothing but the generator of time translations for the vari- 
able r as shown in Ref. 42. In other words, pg d = 
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FIG. 4: Visualizing the modular transformation that maps two copies 
of the solid torus B 2 x S 1 to the pair of interlocking solid tori on the 
right. Expanding the interlocking tori to fill all of space and gluing 
them along their mutual boundary produces S 3 . 



Z _1 exp (—2ttRH t ) with 2ttR playing the role of the inverse 
temperature. If we ignore the conformal factor, as we are al- 
lowed to do for a conformal field theory, then we may compute 
the partition function Z = tr(exp (—2nRH T )) by continuing 
the metric above to imaginary time r = — ite and evaluating 
the partition function of the conformal field theory on the re- 
sulting space. Remarkably, there is a simplification, because 
when we consider the Euclidean version of 10 it can be shown 
to give the round metric on S d+1 . Thus, to obtain the entan- 
glement entropy of a ball of radius R, we need the conformal 
field theory partition function on a sphere of radius R (for de- 
tails see Ref. 42). We also note that in principle we have 
access to the full entanglement Hamiltonian. 

The simplest calculation that illustrates this tool is a com- 
putation of the entanglement entropy of a disk in topological 
liquids in 2 + 1 dimensions. Topological field theories are 
certainly conformal field theories since they do not care about 
the metric of space time at all. To compute Z(S 3 ) we begin 
with the fact that Z(S 2 x S 1 ) = 1. This is the statement that a 
topological phase has a unique ground state on the sphere. The 
spacetime S 2 x S 1 may be cut open along the S 2 to yield two 
copies of the solid torus B 2 xS 1 . Since dB 2 x S 1 = S 1 x S\ 
Z(B 2 x S 1 ) = l^) is a state in the Hilbert space of the torus 
generated by imaginary time evolution. This state is normal- 
ized since 1 = Z(S 2 x S 1 ) = (\f f |\f f ). Instead of gluing 
the tori back together directly, we first make an S modular 
transformation of one of the boundary tori which exchanges 
the two non-contractible surface loops. Now gluing the tori 
together yields S 3 as shown in Fig. 4. The modular transfor- 
mation is implemented using the modular 5-matrix S£, and 
a simple calculation gives Z(S 3 ) = (*|5|*) = Sq. Since 
Sq = 1/2? and using the fact that the ordinary Hamiltonian 
of a topological phase is zero, we conclude that S(B 2 ) = 
\nZ(S 3 ) = — In 2? as originally shown in Refs. 7,8. In fact, 
using the exponentially fast factorization of the density matrix 
for widely separated regions and the topological invariance of 
Aopo' we ma y immediately establish that /?topo = — l n ^ 
where Ng is the number of boundaries of the region of inter- 
est. 



B. General computation 



Using the tools we just introduced, we now give another 
heuristic argument for the conclusion that ftxv = Pxy + 
Pz 2 . Since the XY* critical point is a conformal field the- 
ory, the entanglement entropy of a disk may be obtained from 
Zxy'(S 3 ). We take Z XY (S 3 ) and Z Z2 {S 3 ) as given. To 
obtain XY* from XY ® Z2, we must bind the two sectors 
together by attaching Z2 Wilson lines to the bosons. How- 
ever, in the extreme topological limit we considered above, all 
Wilson lines are tensionless objects in the gauge sector. Thus 
assuming a large gap to visons, the attachment of Wilson lines 
to bosons leaves the boson dynamics essentially unaffected 
as they acquire neither interesting phases nor energetic penal- 
ties. For example, the projected wavef unction we wrote down 
above has the same boson density correlators as the unpro- 
jected wavefunction. Hence we argue that for the purposes of 
universal physics we have Zxy*(S 3 ) ps Zxy(S 3 )Zz 2 (S 3 ). 
In this topological limit we have Z% 2 (S 3 ) = 1/2? =1/2, and 
we recover the fact that (3xy* = Pxy + /3z 2 • 



The advantage of this argument is that it is simple and 
can be easily generalized. Consider any other bosonic criti- 
cal point between an insulating and a broken symmetry phase 
where the bosons are charged under some discrete group G. 
By gauging the group G we can describe deconfined transi- 
tions between exotic insulators and conventional broken sym- 
metry phases. The projection procedure works as above, 
namely, we attach Wilson lines to the bosons according to 
the representation of G they carry. However, even when G is 
non-Abelian, these pure charge Wilson lines (chargons) have 
trivial braiding properties. Furthermore, assuming as before 
that the magnetic excitations (fluxons and flux-charge com- 
posites) are gapped, we can conclude that the boson dynamics 
are not seriously modified by the projection. This allows us to 
conclude that ^ gauge d bosons ~ Z bosons Z G which, via our 
argument above, gives /3 gauged bosons = /3 bosons + Pg- 



We wish to emphasize at this point that the partition func- 
tion argument applies only the special case of the disk. Be- 
cause of the topological nature of /3q, we expect that the disk 
result will hold true for any other simply connected region 
(or multiple such regions). We will simply need to take into 
account the shape dependence of /?g a pless - However, we can- 
not directly probe the topological contribution to more com- 
plicated geometries such as the annulus using this argument. 
Nevertheless, the general considerations laid out above for the 
special case of the XY* critical point also lead us to conjec- 
ture that the topological part of the entanglement entropy for 
the annulus remains at its value inside the topological phase. 
To the extent that the annulus entropy is equal to that of two 
factorized disks, the partition function argument above rigor- 
ously establishes our claim. 
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IV. NON-COMPACT CP 1 MODEL 

A. Description of transition 

Here we consider certain square lattice spin- 1/2 quantum 
antiferromagnets realizing a direct continuous transition be- 
tween a Neel ordered state and a valence bond solid (VBS) 
state. Since these two phases break different symmetries 
present in the microscopic Hamiltonian, a continuous transi- 
tion between them is non-generic in the Landau theory. Un- 
less we fine tuned the model, we would expect either a first or- 
der transition or coexistence. Thus the critical point, if there 
is one, cannot be the conventional 0(3) critical point. The 
mechanism for the phase transition is actually the condensa- 
tion of spinons, bosons with fractional spin which are confined 
in the Neel and VBS phases, but which are liberated at the 
critical point. The low energy description of the transition is 
given in terms of the so-called non-compact CP 1 model. This 
model consists of bosonic spinons carrying fractional charge 
coupled to a fluctuating gauge field. Right at the critical point, 
the instantons which typically confine a U(l) gauge theory in 
2 + 1 dimensions become irrelevant. Lattice models that real- 
ize this transition differ from the pure field theory by irrelevant 
operators that flow to zero in the low energy limit. 

B. Irrelevant operators 

Having described the basic physics of the Neel to VBS tran- 
sition and its potential description in terms of the non-compact 
CP 1 model, we would like to answer the following question. 
Given that monopoles are irrelevant at the critical point, so 
that the gauge field is effectively non-compact, will a field the- 
ory computation, say using continuum QED, give the correct 
universal features of the entanglement entropy? For example, 
will (3 as defined by the lattice model at its critical point be the 
same as (3 computed using the field theory. We argue that this 
is so since these two models differ only by irrelevant operators 
at the critical point. In particular, the microscopically allowed 
instanton operators are all irrelevant at the critical point. Of 
course, away from the critical point there will be a compli- 
cated crossover description, for example, the Neel phase will 
have its own value of (3 due to gapless spin waves. 

Another way to see this is as follows. We can always for- 
mulate the transition in terms of the VBS order parameter so 
long as we include the appropriate WZW terms responsible 
for adding quantum numbers to the vortices of the VBS order 
parameter (see Refs. 43,44). The presence of monopoles is 
related to the four fold anistropy of the VBS order parame- 
ter. Since this anistropy is irrelevant at the critical point, we 
may as well study the model with full SO (2) symmetry in the 
presence of the WZW term. This model may again be shown 
to be equivalent to the non-compact gauge theory via a duality 
transformation. 

Note that this is not conflict with the comment made in the 
introduction about the importance of remembering the com- 
pactness of order parameter fields. For example, we should 
not expect f3 in the superfluid phase to agree with (3 of a free 



boson, as one is compact and the other is not. The essential 
point is that the compact model may be viewed as a gauged 
version of the non-compact model, and this gauging proce- 
dure means that many operators that are apparently sensible in 
the low energy theory are in fact not allowed. Related to this, 
there must also be new operators, vortex operators, that corre- 
spond to the topological defects of the model. In the case of 
the gauge description relevant for the Neel to VBS transition, 
the low energy theory already contains all the fields needed to 
insure that no unphysical operators can be constructed. Our 
claim is thus that the compactness of a U(l) gauge field can 
be forgotten when monopoles are irrelevant even though the 
compactness of an order parameter cannot. 



C. Renormalization group flow 

Besides arguing for the insensitivity of the universal parts 
of the entanglement entropy on irrelevant operators, we would 
also like to comment on the extent to which "exotic" crit- 
ical points may be regarded as more entangled that simple 
symmetry breaking critical points. By analogy with topolog- 
ical phases in two dimensions, we might suggest that a larger 
(3 corresponds to more entangled critical point. Certainly it 
seems that phases with greater long range entanglement and 
more complicated topological order generally have a larger to- 
tal quantum dimension. However, we must also be careful to 
note that we can typically increase (3 by simply tensoring to- 
gether many weakly interacting degrees of freedom. The most 
sensible comparisons of f3 between critical points would thus 
be made when those critical points describe transitions out of 
the same phase, so that we in some sense fix the number of 
degrees of freedom. 

For example, we could consider two transitions out of the 
Neel ordered state. The first is a conventional symmetry 
breaking transition while the second is the deconfined criti- 
cal point we have been considering. If j3 were larger for the 
deconfined critical point, then we would be justified in saying 
that it was more entangled than the conventional 0(3) criti- 
cal point. Note that we might in general have to be careful 
about choosing the right region when trying to establish such 
a claim. Below we compare f3s for disks. 

In fact, a conjecture is available that would imply the decon- 
fined critical point is more entangled, suitably defined, than 
the conventional 0(3) critical point. In Ref. 32 it is proposed 
that a certain subleading term in the entanglement entropy of 
half of a sphere is monotonic under RG flow. This claim is 
supported with a variety of holographic calculations, and the 
claim in 1 + 1 and 3 + 1 dimensions appears to be established 
thanks to c-theorems in those dimensions, see Ref. 45. As- 
suming this claim is true in 2 + 1 dimensions, we can say that 
the deconfined critical point is more entangled than the 0(3) 
critical point because it can flow to the usual 0(3) point un- 
der relevant deformations 2 . To be precise, if the full system is 



2 For example, the single instanton operator is relevant at the non-compact 
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tuned to criticality on a sphere, then the entanglement entropy 
of half the sphere contains a subleading term that decreases 
under RG flow from the deconfined critical point to the 0(3) 
critical point. We must be careful in extending these claims 
too far, however, since the results of Refs. 46-48 show that (3 
can be negative for some Lifshitz critical points. It is not clear 
how to reconcile this behavior with a general entanglement 
c-theorem. 



V. EXTRACTING UNIVERSAL TERMS AND NUMERICS 
A. General considerations 

We now return to the somewhat subtle issue of properly 
defining and extracting the universal coefficient j3. In practice, 
the simplest method is simply to compute S(L) for a variety of 
region sizes and fit to S(L) = aL — (3+0(l/L). This method 
seems to work well enough, but it is dangerous because of 
the possibility of mixing some of the non-universal compo- 
nent in f3. This is especially problematic when the region is 
question has corners, as it often does in a lattice formulation. 
Corners are expected to give a contribution like log (L/e) and 
thus changes in e, the UV regulator, can look like changes in 
the "universal" parameter f3. It is therefore important to have 
a well defined way, at least in principle, to extract j3. 

For topological phases in two dimensions the constant /3 
is related to the total quantum dimension of the topological 
phase. It can be measured independent of any non-universal 
components of the entanglement entropy using a particular set 
of subtractions that we have already discussed. This kind of 
procedure is familiar from the theory of mutual information 
X(A, B) = S(A) + S(B) - S{AB) where the subtractions 
can be motivated in the continuum by the desire to remove the 
non-universal boundary law terms (singularities still arise if A 
and B touch, see Ref. 28). The finite correlation length in a 
gapped phase both implies that distant parts of the boundary 
may be treated separately and that corners do not give special 
contributions. Both of these statements may fail in a gapless 
phase. 

As already mentioned, corners are known to give a logarith- 
mic correction to the entanglement entropy. This correction 
is analagous to the singularities that appear in Wilson lines 
with sharp corners in gauge theory. From the point of view of 
the renormalization group argument sketched in the introduc- 
tion, the logarithm arises because the corner contributes order 
one entanglement at every RG step. Thus if we integrate the 
corner contribution from e to L we find a contribution to the 
entanglement entropy going like J dr/r ~ log (L/e). How- 
ever, one can check that the substraction used to extract j3 for 
topological phases also correctly removes any potential cor- 
ner terms. This works in part because the coefficient of the 
logarithmic term is an even function of the angle measuring 
the deviation from straightness i.e. a 7r/2 corner gives the 



CP 1 critical point with the RG flow terminating at the 0(3) critical point. 



same contribution as a 3ir/2 corner. On the practical side, the 
contributions from corners can often be quite small. It can 
be further mitigated by considering very slight corners, how- 
ever, such an arrangement often requires a very complicated 
region shape with many small corners. It is possible that the 
best proposal is to simply use a "rough" circle on the lattice 
with many small lattice scale corners trusting the RG flow to 
quickly smooth out these irregularities. 

The other potential concern is the slow (power law) fall-off 
of correlations in a gapless phase. For example, the subtrac- 
tion to extract j3 relies on the assumption that the subleading 
term in the entanglement entropy of an annulus is — 2P^ s u. 
However, we can already see a problem with this assumption, 
namely that f3 can have non-trivial shape dependence, so we 
must ask if the annulus is — 2/3 ( }j s j c or something else entirely. 
It is clear that for widely separated regions the entanglement 
entropies will add, however the annulus geometry is not ex- 
actly of this type. We may gain some intuition for this situa- 
tion by appealing to the twist operator formalism for entangle- 
ment entropy developed in Ref. 28. That formalism suggests 
the following test: if the correlation functions of twist fields 
decay quickly enough then the inner and outer ring of the an- 
nulus may be treated separately (up to power law decaying 
corrections). To estimate their decay, we use the fact that at 
large separation the mutual information (which involves in- 
tegrating the twist field correlator over distant region bound- 
aries) bounds the connected two point function of local oper- 
ators. If Ao is the smallest scaling dimension among all local 
gauge-invariant operators, then we expect T(x) ~ 1 /x 4A " at 
large separation x. The unitarity bound in 2 + 1 dimensional 
conformal field theories requires that Ao > 1/2 for scalar op- 
erators. Practically speaking, the corrections could be quite 
large for a finite annulus, so one would have to study a range 
of annuli sizes to extract the constant term from finite size 
scaling. 

All these considerations lead to a much more complicated 
situation in the case of gapless phases. However, we are still 
interested in studying the entanglement properties of critical 
points in terms of (3 and corner contributions. The simplest 
protocol is thus to use smooth simply connected regions to 
extract (3 from finite size scaling. One could further study the 
shape dependence of j3 as a diagnostic of our proposed decom- 
position into gapless and topological contributions (although 
this decomposition may not hold in more complicated theories 
where the gapless states are not pure charge). Unfortunately, 
this may be nearly impossible in a lattice setup where corners 
are almost inevitable. We could still isolate the shape depen- 
dent part of (3 by studying, for example, Ssquare — ^rectangle- 
This subtraction will approximately remove corner terms, and 
if we keep the perimeter of the rectangle equal to that of the 
square, the subtraction will depend only on the aspect ratio 
of the rectangle. However, this is still not fully satisfactory 
since the topological part of j3 will simply be subtracted away. 
Another possibility is to consider the system on a torus and 
study regions that wrap around the loops of the torus. In this 
way it is possible to have a region on a lattice without corners, 
but at the cost of exposing the entanglement entropy to addi- 
tional topological subtleties. However, at least in the case of 
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gapped topological phases, these subtleties can be understood 
and taken into account. 

Let us summarize the situation as we see it. If we do use the 
usual topological protocol for gapless systems, then we will 
obtain the topological part plus a sum of terms that depend on 
region shape. It is possible that these shape dependent terms 
can be calculated or fit numerically and subtracted away. The 
result would be the topological part of the entropy, however, 
we stress that this would be a quite complicated procedure in 
general. Corners are present and are hard to avoid on the lat- 
tice. Regions wrapping around the entire compact space can 
avoid corners but come with additional complications. There 
will be power law corrections to the various asymptotic state- 
ments we have made in this paper, as well as the aforemen- 
tioned shape dependent universal terms coming from the gap- 
less modes. Results from field theory often place special em- 
phasis on the properties of simple spaces like the sphere and 
simple regions like disks, regions that are hard to replicate on 
the lattice. Finally, we remind the reader again that finite size 
corrections can have a very profound impact on the entangle- 
ment entropy e.g. freezing out gapless degrees of freedom. 
Despite all this complexity, as we have illustrated in this pa- 
per, it is possible to gain control of the universal parts of the 
entanglement entropy in a wide variety of geometries. Thus 
we believe we have the necessary flexibility to untangle the 
complexity of universal gapless entanglement. 

B. Numerical issues 

In the numerical investigation of the XY* transition 
18, there were some puzzles associated with the observed 
entropies. Those authors considered a hard core boson 
model with nearest neighbor hopping on the kagome lat- 
tice with an extra constraint term associated with hexagons. 
The Hamiltonian is H = —t^ <rrl> b+b r i + b^,b r + 
V ShexEiehex n i) 2 - This model has a superfluid to Z2 in- 
sulator transition with increasing V located at V/t ps 7. In 
the absence of V the model is in a superfluid phase where the 
elementary excitations are phonons. The speed of these ex- 
citations will be set by t in units of the lattice spacing since 
there is no other scale. To probe the contributions to f3 from 
critical degrees of freedom we must consider a system of large 
enough size and low enough temperature so that these low en- 
ergy degrees of freedom can fluctuate. At the critical point, 
similar considerations apply, and although the velocity can be 
renormalized, we still expect relativistic low energy modes 
with speed ~ V ~ t. Given the lattice sizes, region sizes, and 
temperatures considered in Ref. 18, it is possible that some of 
these gapless modes were frozen. Thus we cannot be sure that 
the shape dependent part of the entropy was being effectively 
sampled. 

We have also predicted a topological part at the critical 
point (as measured by the usual protocol) given by the full 



—2 In 2 of the insulating phase. This is so despite the pres- 
ence of fluctuating Z2 charges as we argued above. However, 
it is also possible, depending on the details of the temperature 
and the energy scale of the charge motion, that these fluctu- 
ating charges are moving incoherently, a situation that then 
gives a topological part of — In 2. If we are really probing the 
ground state properties, then we believe the topological part 
should be —2 In 2. Even if this number is not observed ex- 
actly, the fact that topological contributions to the entropy are 
not continuously variable means we can still be confident we 
have a Z2 deconfined critical point or phase. 

The results of Ref. 18 are quite interesting, and more work 
is needed on both sides to achieve a completely satisfactory 
comparison between theory and experiment. We hope that the 
successes of Ref. 18 lead to other calculations of the entan- 
glement entropy at quantum critical points, conventional or 
otherwise. 



VI. CONCLUSIONS 

We have studied the entanglement structure of a variety of 
deconfined quantum critical points focusing on 2 + 1 dimen- 
sions. We have computed the universal terms in the entangle- 
ment entropy in terms of the same data at conventional criti- 
cal points and in topological phases. We have also argued for 
the insensitivity of the universal terms to irrelevant operators 
in the low energy description. Based on our explicit calcu- 
lations and on a conjectured c-theorem for entanglement en- 
tropy, we can confidently assert that deconfined critical points 
are more entangled, in the sense discussed above, than cor- 
responding conventional critical points. Finally, we carefully 
analyzed the structure of entanglement for more complicated 
region geometries and discussed in detail various protocols for 
extracting universal terms. 

Given the power of the partition function argument for com- 
puting the entanglement entropy of balls in a CFT, there are 
many interesting systems that could be attacked using this 
technology. It would also be very interesting to investigate 
further the full entanglement Hamiltonian in a CFT, which is 
known for some special regions. Despite the boundary law 
for entanglement entropy, the entanglement Hamiltonian in a 
CFT is not well localized to the boundary, so it would be in- 
teresting to understand more precisely how local the entangle- 
ment Hamiltonian is. There is also considerable opportunity 
to carry out further numerical studies of entanglement in de- 
confined critical points (and, in fact, in conventional critical 
points). Thanks to powerful new perspectives and tools, there 
is a growing body of analytical results in higher dimensions 
that can be compared with numerical calculations. 
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